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Abstract 

We consider an exactly tractable model of the Kramers type for the voltage-dependent 
gating dynamics of single ion channels. It is assumed that the gating dynamics is 
caused by the thermally activated transitions in a bistable potential. Moreover, the 
closed state of the channel is highly degenerate and embraces the whole manifold 
of closed substates. Opening of the ion channel is energetically prohibited from 
most of the closed substates and requires a special conformation where the voltage 
sensor can move along an activation pathway and trigger the transition into the 
open conformation. When the corresponding activation barrier towards the chan- 
nel's opening is removed by the applied voltage, the statistics of non-conducting 
time intervals become strongly influenced by the conformational diffusion. For the 
corresponding supra-threshold voltages, our model explains the origin of the power 
law distribution of the closed time intervals. The exponential-linear dependence of 
the opening rate on voltage, often used as an experimental fit, is also reproduced 
by our model. 
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1 Introduction 



Ion channels present complex protein structures embedded across the biolog- 
ical cell membranes thereby forming the ion-conducting, selective nanopores 
[1]. The conformational dynamics of these special proteins, which is known 
under the label of gating dynamics, results in the spontaneous openings and 
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closures of ion channels [1]. In the simplest case, this gating dynamics has 
the "all-or-none" character and can be symbolized by the following kinetic 
scheme, 



C 



O . 



(1) 



Accordingly, the conductance of ion channel fluctuates stochastically between 
some finite value and nearly zero. In the case of voltage-dependent ion chan- 
nels, like a delayed rectifier potassium channel, or a Shaker K"*" channel, the 
opening rate, koiV), and the closing rate, kciV), are both dependent on the 
applied voltage V. 

The invention of the patch clamp technique [2] enables one to observe the cur- 
rent and the conductance fluctuations in real time with a resolution on the level 

of single ion channels. Moreover, the study of the statistics of dwelling time- 
intervals becomes feasible. As a matter of fact, the patch clamp experiments 
with single ion channels pioneered the whole area of single molecular research. 
The gating dynamics of an ion channel with one closed and one open states, 
as seen visually in the patch clamp recordings, can be characterized by the 
probability densities of closed, /c(t), and open, /o(t) dwelling time- intervals. 
The experimental studies reveal that in many ion channels the distribution of 
open dwelling times is approximately exponential, fo{t) = kcexp{—kct). How- 
ever, the distribution of closed time intervals fc{t) frequently involves the sum 
of many exponentials, Aiexp(— Ajt), i.e.. 



with weights Cj, obeying Ci — 1. The rationale behind this fitting proce- 
dure is the assumption that the closed state consists of N kinetically distinct 
discrete substates separated by large potential barriers. The gating dynam- 
ics is then captured by a discrete Markovian kinetic scheme with the rate 
constants determined by the Arrhenius law. 

Such an approach presents nowadays the standard in the neurophysiology [3]. 
An essential drawback of it is, however, that the number of closed substates 
needed for agreement with the experimental data can depend on the range of 
applied voltages and temperatures used in the experiment. For example, the 
experimental gating dynamics of a Shaker potassium channel has been reason- 
ably described by a sequential 8-state Markovian scheme with 7 closed states 
for a temperature about T = 20 °C [4]. However, when the same research 
group extended later on their model to embrace the small range of tempera- 
tures 10 — 20 "C, three additional closed states have been introduced [5]. This 



N 



fc{t) = X] Ci\ exp(-Aii), 
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"proliferation" of the number of discrete states, required for the consistent de- 
scription of experimental data, is rarely addressed within the approach based 
on a discrete Markovian description. 

Moreover, it may happen that for some channels the closed time distribution 
can be more conveniently fitted by a stretched exponential [6], or by a power 
law dependence [7,8,9] with a few parameters only. This observation gave 
rise to the so-called fractal models of the ion channel gating dynamics such 
as put forward, by example, by Liebovitch et al. [6]. The diffusion models 
introduced by Millhauscr et al. [7], Lauger [10] and Condat and Jackie [11] 
are intrinsically based on the concept of dynamical conformational substates 
in proteins, an idea which has been pioneered by Frauenfelder et al. [12]. The 
diffusion models serve as a microscopic justification for the fractal modeling. 
On the other hand, the discrete diffusion models yet present Markovian models 
with a large number of states. The non-Markovian, fractal behavior emerges 
from the projection of a full Markovian dynamics onto the subspace of two 
observable states as symbolized by the kinetic scheme (1). 

Alternatively, the discrete diffusion can be replaced by the continuous diffu- 
sion on a potential landscape. Then, the distinct minima of this landscape, 
separated by substantial potential barriers, correspond to the discrete states 
in the mainstream approach. However, it may happen that, depending on the 
applied voltage, some of the distinct barriers disappear. Then, the new fea- 
tures of gating dynamics come into the play, which are inevitably missed in 
the discrete modeling. It seems therefore that a compromise between the dis- 
crete Markovian schemes and a continuous conformational diffusion approach 
is called for [13,14]. Especially, the continuous diffusion models, if set up suf- 
ficiently simple to allow an analytic treatment, are capable to provide a new 
insight into the problem of ion channel gating. In the present work, we refine 
and justify further the approach put forward in Ref. [14]. 



2 Model of gating dynamics 

To start, let us consider the model depicted in Fig. 1. Its central element 
is the volt age- dependent transitions in a bistable, piecewise linear potential 
along the x-axis, see l.h.s. in Fig. 1. This bistable potential corresponds to the 
motion of the so-called voltage sensor. The voltage sensor presents a functional 
part of the ion channel machinery which changes its position in response to 
the changes in the transmembrane potential. In accord with the current view 
[1], the voltage sensor is formed by the system of four positively charged S4 
a-helices which carry a total gating charge g ~ 10 — 13 Cq, where cq is the 
positively-valued, elementary charge. When a; = 0, the voltage sensor is in its 
resting position; when x — Xq, the voltage sensor is fully activated and the 
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Fig. 1. Sketch of the studied model. The left part of figure involves the 
two-dimensional conformational space for the gating dynamics. The x-axis corre- 
sponds to the voltage sensor degrees of freedom and the y-axis to the conformational 
diffusion. Note that only the y = cross-section, U{x,0), of the two-dimensional 
conformational potential U{x, y) possesses a bistable character. For y > 5 exceeding 
some small 5-vicinity of y = 0, the y = a = constant cross-sections, U{x,a), are 
essentially monostable, exhibiting sharp minima at a; = 0. Moreover, the potential 
U{0, y) is flat in the y-dircction. These features allow one to simplify the picture 
of a two-dimensional reaction state space to the case of one-dimensional reaction 
coordinate (see the right part in the figure) as described in the text. 

complete gating charge is moved across the membrane. Moreover, the drop 
V of the transmembrane electric potential is assumed to occur linearly on 
the characteristic length xq which corresponds to the effective width of the 
cellular membrane within the region of the gating charge translocation. Then, 
the energy barriers, Uo{V) and Uc{V), follow as 



Uo{V) = qr{V,~V), (3) 
C/e(l^) = C/f -g(l-r)(K-^), (4) 

where r := Xb/xq. Furthermore, Vc in (3) corresponds to the threshold value 
of the voltage V when the activation barrier towards the channel opening 
vanishes and U^^^ in (4) is the corresponding height of the activation barrier 
towards the channel closing. Note that for V > Uo{V) is negative, but 
Uc{V) has always a large positive value, since one assumes that U^^^ is large 
and 1 — r << 1. These assumptions correspond to the experimental observa- 
tions that the closing rate has a very strong dependence on the temperature 
and it is less voltage-sensitive than the opening rate [1]. Furthermore, one as- 
sumes that the voltage sensor located at the position x > xb {xb corresponds 
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to the top of the activation barrier for the sub-threshold vohages, V < Vc) 
triggers the conformation change in the activation gate. This latter confor- 
mational change then leads to the channel opening. It is worth to notice that 
in accord with the reasoning in [15,16], the activation gate is likely formed 
by the bundle of the intracellular ends of inner, pore- lining S6 a-helices [1]. 
The motion of the voltage sensor creates an indirect mechanical torque on the 
activation gate. This torque is mediated by other structural elements of the 
channel protein. For this reason, the link between the voltage sensor and the 
activation gate may be flexible and this circumstance can introduce a kind of 
mechanical hysteresis. Namely, the closure of the activation gate does not oc- 
cur immediately, when the voltage sensor crosses the border at a; = a;^ in the 
back direction, but rather requires that (is most probable when) the voltage 
sensor is fully returned to its resting position at a; = 0. For this reason, by 
calculating the dwelling time distribution fc{t) we assume the initial condition 
for the probability density P{x,t) in the form P(a;, 0) = 5{x), i.e., each and 
every closing time-interval starts when the voltage sensor has fully returned. 
This inevitably presents an approximation to the reality. In accord with the 
general activation theory of Kramers [17], we also have to impose an absorbing 
boundary at x > xb- The most natural choice for this is x = Xq. However, in 
view of our model assumption Uf''' » ksT one can safely move the absorb- 
ing boundary to x = xb, i.e., P{xB,t) = at all times t > 0. The reason 
for this is that the final transition xb ^ xq leading to the opening is very 
quick. The same is assumed in many discrete Markovian schemes [1]. Given 
this assumption, the duration of this final step is negligible in calculating fc{t). 
Furthermore, in order to find the distribution of open times, fo{t), one should 
put the refiecting boundary at x = xq and the sink (absorbing boundary) at 
X — 0, in accord with the above discussion. 

The next important element of our model is the assumption that the voltage 
sensor can operate only if the protein acquires a special pre-open conforma- 
tion. Otherwise, its motion is energetically prohibited due to the steric factors. 
To account for this type of behavior we introduce the additional, "perpendic- 
ular" conformational reaction coordinate y, see l.h.s. in Fig. 1. It describes 
the conformational diffusion. One assumes that the corresponding dynamics 
is only weakly voltage-dependent and we shall neglect the accompanying re- 
distribution of charges in the protein molecule. Actually, this conformational 
diffusion should occur in a hierarchically arranged, ultrametric conformational 
space [12,18,19]. In a simplified manner, one can model this hierarchical fea- 
tures by a random walk with the mean-squared displacement {Ay^{t)) oc t^. 
This latter relation holds, strictly speaking, if the diffusion is unbounded in 
space. The bounded character, however, is important for the following. We 
restrict ourselves in this paper to the case of normal diffusion, 7 = 1, as the 
simplest possibility. Finally, within this outlined picture one can "convert" 
the dynamics along the |/-axis into the dynamics along the x-axis, by use of 
the extra part [— L, 0] there, see r.h.s in Fig. 1, whereby introducing the aux- 
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iliary diffusion length L and the diffusion coefficient D = ksT/ri. Here, 7] is 
the viscous friction coefficient characterizing the conformational dynamics. In 
order to account for the bounded character of the conformational diffusion, 
we impose the reflecting boundary condition ^^^^\x=-l — ^ sX x — —L for 
all t. Note that 77 should generally depend exponentially on the temperature, 
i.e. 

77 = r]oexp{e/kBT), 

where e is a characteristic activation energy for transitions between conforma- 
tional microstates. It can be in the range of several fcsTroom- The parameters 
D and L are conveniently combined into the single parameter - the confor- 
mational diffusion time, i.e.. 

To = L^D oc eM^/kBT)/kBT, (5) 
which constitutes an important parameter of the theory. 



3 Theory of ion channel gating 



In order to find fc{t) one has to solve the Smoluchowski dynamics 



where f3 = l/lksT), supplemented by the initial and the boundary condi- 
tions of reflection and absorption discussed above. The closed residence time 
distribution then follows as 

m - (7) 

where ^c{t) — /ff P{x, t)dx is the survival probabihty in the closed state. 

By use of the standard Laplace transform method we arrive at the following 
exact solution for the Laplace-transformed distribution of closed times /c(s), 
see in Ref. [14]: 

iM ^ (8) 

where 
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A{s)^eM-(3Uo{V)l2)y[^uf^^)^^Ae^s (9) 
B{s) = ^ ^W^,{V) + ^erDS cosh (^^l3Wl{V)+AeTDs/2^ (10) 

+ (2eV^tanh ^ - (3Uo{V)j sinh <^^J (3Wl{V) + AeTDsl2}j , 

where the parameter C, is given hj C, '■= Xb/L. The expUcit resuh in Eqs. (8)- 
(10) allows one to find all moments of the closed residence time distribution. 
In particular, the mean closed residence time, 

(T,)- / tMt)dt = lmi[l-Ms)]/s, 

Jo s^O 

emerges as 



^ , Wl-)(e/^^°(^)-l-0+e(e^^°(^)-l) 

The effective opening rate can be defined as ko{V) :— l/{Tc{V)). Let us 
consider the limiting case ^ = xb/L -C 1. In the language of discrete sub- 
states, this limit is tantamount to the assumption that the number of quasi- 
degenerate conformational substates, which correspond to the resting position 
of the voltage sensor, largely exceeds that of the voltage-sensor. Under this 
assumption, we obtain in leading order of ^ 

Note that the functional form in Eq. (12) is nothing but the familiar exponential- 
linear dependence used as a phcnomcnological experimental fit in the cel- 
ebrated paper by Hodgkin and Huxley [20] used to describe the voltage- 
dependence of the opening rate of a single gate in the potassium channels. 
This form is commonly used to parameterize the opening rate of the potas- 
sium channels, see, e. g., in [21]. Our model provides a detailed justification 
for this experimental result. Its remarkable feature is that the dependence 
of the rate on voltage is exponential for V < Vc, when the energy barrier 
towards activation of the voltage sensor is essential, ksT < Uo{V). This expo- 
nential voltage-dependence implies in virtue of (5) an exponential dependence 
on temperature as well, i.e., 

koOiexp{-[e + rq{Vc-V)]/kBT} . (13) 

This exponential temperature dependence has two sources: a first one is due 
to the activation barrier of the voltage sensor Uo{V), while a second one is 
due to the activation barrier e between diffusional microstates, which we have 
assumed for reasons of simphcity to be of equal height. The barrier Uo{V) 
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Fig. 2. Closed time probability density from our model, Eq. (7), (solid line) and 

6-exponential fitting procedure of the experimental data from Ref. [8] (dotted line). 
The following set of parameters is used in our calculations: td = 2.9 sec, $, = 0.01 
and f5Uo{V) = -1.653. 

can acquire large values. For example, assuming typical values Vc = —40 
meV, rq — 10 Cq and a room temperature kBTj-oom — 25 meV one obtains 
UoiYr) — 20 ksT for the cell resting potential = —90 meV. Furthermore, 
when the activation barrier UoiV) vanishes for V > Vc-, 'ii follows from (12) that 
the rate koiV) exhibits the linear dependence on voltage, i.e., ko{V) oc {V —Vc). 
In this case, its temperature dependence is distinctly reduced and becomes 
mainly determined by the activation barrier e of the conformational diffusion. 
The latter one can assume a few kBTroom only. The very different temperature 
dependences of the opening rate for V ^ Vc and for y > present an 
interesting feature of our model which calls for an experimental verification. 

It is worthwhile to notice that the discussed crossover from an exponential 
to linear voltage dependence of the opening rate is qualitatively preserved for 
any value of C,, including the extreme case ^ ^ 1. In this case L <^ Xb and 
the conformational diffusion does not play a dominant role. The discussed 
feature is caused by the vanishing of the activation barrier Uo{V) for V > Vc- 
However, the distribution of closed time-intervals fc{t) depends qualitatively 
on ^. Namely, for ^ <^ 1 and V >Vc\t displays a power law regime, i.e. 



2{nTDy/^kot^/^' 



(14) 



for the intermediate time scale {Tc)'^/td <^ t <^ td. In Fig. 2, we illustrate this 
prominent feature for the following model parameters: td — 2.9 sec, ^ = 0.01 
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Fig. 3. Logarithmic derivative of the closed time probability density fc{t)- The solid 
line and the dotted line correspond to those in Fig. 2; the long-dashed level line 
corresponds to strict power law dependence fc{t) oc The oscillating behavior 

reveals the hierarchical character of the conformational dynamics. 

and f3Uo{V) = —1.653. The solid line denotes the closed time probability den- 
sity obtained from our model, Eqs. (8)- (10), by a numerical inversion of the 
Laplace transform fc{s). The short-dashed line presents the fitting of the ex- 
perimental data for the delayed rectifier K"*" channel from a neuroblastoma x 
glioma cell by use of Eq. (2) with 6 exponentials. This fitting procedure is taken 
from Ref. [8] (see, Table 4 therein) and implicitly corresponds to a discrete 
Markovian scheme with 6 closed substates. Both results describe well the in- 
termediate power law trend and the exponential tail of the closed time-interval 
distribution. Nevertheless, some small distinctions can be distinguished. 

The plot of the logarithmic derivative, dlog fc{t)/d log t, versus the logarithm 
of time t in Fig. 3 renders these distinctions much more visible. The advantage 
of such a representation as in Fig. 3 is - in accordance with the reasoning in 
Ref. [19] - that the hierarchical, tree-like relaxation dynamics reveals itself by 
logarithmic oscillations around the level line that corresponds to the power 
law trend. Remarkably enough, the exponential fit from Ref. [8] does exhibit 
such logarithmic oscillations, cf. Fig. 3. Thus, this result seems to support the 
hierarchical character of the conformational dynamics. Our simplified model 
does not distinctly display these fine features as these are rooted in the dis- 
crete nature of hierarchical states. Nevertheless, the power law trend, which 
reveals the presence and the importance of the conformational dynamics, is 
reproduced by our model. Moreover, its appealing feature is that it requires 
only few parameters which possess a clear physical meaning. The particular 
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value of the power law exponent —3/2 corresponds to the conformational dy- 
namics modeled as a bounded normal diffusion. Other power law exponents, 
also seen experimentally [9], require a generalization of our model to the case 
of anomalous diffusion. Such corresponding work is presently in progress. 



4 Summciry and conclusions 

We herewith have presented a simple model of the complex gating dynamics 
in voltage-dependent potassium ion channels. It is based on the concept of 
conformational diffusion. In particular, we assumed that the conformational 
change leading to the opening of ion channel is triggered by the voltage sensor 
which, in its turn, is permitted only when the channel protein acquires a special 
pre-open configuration. When the ion channel is closed, it exhibits an inter- 
nal, conformational diffusion over the manifold of conformational substates 
which do not possess a sensitive voltage dependence. In a simplified manner, 
this conformational diffusion has been mathematically modeled by bounded 
normal diffusion. Moreover, it has been assumed that the open conformation 
of the channel is more ordered, with less conformational substates. Then, the 
conformational diffusion does not play an essential role. In the language of 
statistical thermodynamics this means that the ion channel upon opening un- 
dergoes a kind of ordering transition into a state with lower configurational 
entropy. 

We should also stress here that our simple model is aimed not to replace the 
standard discrete Markovian modeling [3], but rather to complement these 
efforts by highlighting some basic physical features which otherwise become 
blurred with the standard approach. In particular, it has been shown that the 
transition from an exponential to linear voltage dependence of the opening rate 
occurs when the activation barrier for the voltage sensor towards the channel's 
opening vanishes due to the applied transmembrane voltage. Moreover, if the 
conformational diffusion time Tu exceeds largely the mean duration of closed 
time intervals (Tc), the closed time distribution exhibits a power law feature 
on the intermediate time scale {T^'^ /td <^ t <^ td- This power law changes 
over into an exponential tail for times t > td- These features are seemingly 
consistent with the experimental observations for some K+-channels. The true 
physical benchmark of our model is, however, the prediction that the opening 
rate will become much less temperature-dependent for supra-threshold volt- 
ages. The weak temperature dependence in this latter regime should correlate 
with a weak voltage dependence. This distinct prediction calls for an experi- 
mental verification, which hopefully will support our reasoning for the gating 
dynamics in ion channels. 
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